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High-energy Li-Be-B nuclei in cosmic rays are being measured with unprecedent accuracy by the 
AMS experiment. These data bring valuable information to the cosmic ray propagation physics. 
In particular, combined measurements of B/C and Be/B ratios may allow to break the parameter 
degeneracy between the cosmic-ray diffusion coefficient and the size of the propagation region, 
which is crucial for dark matter searches. The parameter determination relies in the calculation 
of the Be and B production from collisions of heavier nuclei with the gas. Using the available 
cross-section data, I present for the first time an evaluation of the nuclear uncertainties and their 
impact in constraining the propagation models. I found that the AMS experiment can provide tight 
constraints on the transport parameters allowing to resolutely break the degeneracy, while nuclear 
uncertainties in the models are found to be a major limiting factor. Once these uncertainties are 
accounted, the degeneracy remains poorly resolved. In particular, the Be/B ratio at ~1- 10 GeV/n 
is found not to bring valuable information for the parameter extraction. On the other hand, precise 
Be/B data at higher energy may be useful to test the nuclear physics inputs of the models. 

PACS numbers: 98.70.Sa, 96.50.sb, 25.40.Sc, 95.35.-l-d 


I. INTRODUCTION 

Understanding the cosmic ray (CR) transport pro¬ 
cesses in the Galaxy is a major subject in modern as¬ 
trophysics. The CR transport is studied using data on 
secondary nuclei {e.g., Li-Be-B), which are created by 
fragmentation of heavier elements, and primary nuclei 
{e.g., H, He, C-N-0), which are produced and acceler¬ 
ated in Galactic sources. Secondary-to-primary ratios, 
and in particular the B/C ratio, are used to constrain 
the Galactic diffusion coefficient, D, and the half vertical 
extent of the propagation region, L [1-3]. For a rigidity- 
dependent diffusion coefficient D k. DqTZ^ , the B/C ratio 
fixes both S and the Dq/L ratio. The degeneracy be¬ 
tween Dq and L may be lifted using data with unstable 
isotopes, such as the ^°Be/®Be isotopic ratio or the Be/B 
elemental ratio [4, 5]. Hence, combined measurements of 
B/C and Be/B ratios may provide the determination of 
the basic CR transport parameters. Understanding CR 
transport is crucial to reliably predict the secondary pro¬ 
duction of antimatter and to set stringent limits on dark 
matter annihilation signals. The parameter L is of great 
importance for assessing the dark matter signal. 

The spectra of B and Be nuclei are now being measured 
by the Alpha Magnetic Spectrometer (AMS) experiment 
in the International Space Station (ISS). Recent mea¬ 
surements of CR protons [6] and preliminary results on 
light nuclei [7, 8] show that AMS is probing the GeV- 
TeV energy region to a ^ % level of accuracy. With this 
standard of precision, it is now timely to assess the theo¬ 
retical uncertainties of the model predictions. In particu¬ 
lar, calculations of Be-B production rates rely on several 
cross-section (XS) estimates. Propagation models make 
use of semi-empirical XS formulae calibrated to acceler¬ 
ator data. Thus, the accuracy of the inferred transport 


parameters is directly linked to the quality of the avail¬ 
able measurements on nuclear fragmentation. 

The aim of this paper is to estimate the nuclear uncer¬ 
tainties in CR propagation and to investigate how they 
affect the determination of the CR transport parameters. 
In particular, I will focus on the anticipated AMS data 
on the B/C and Be/B ratios and their connection with 
the Dcj/L degeneracy problem. For this purpose, I have 
gathered all available XS data for Be and B production 
from B-C-N-0 collisions off hydrogen target. These data 
have been used to constrain the XS parameterizations 
in order to estimate their uncertainties. The resulting 
XS errors have been therefore converted into theoretical 
uncertainties for the model predictions and, finally, into 
uncertainties on the transport parameters that can po¬ 
tentially be inferred by AMS. 


II. CALCULATIONS 

The CR propagation model — This work relies on 
the diffusive-reacceleration model implemented under the 
code GALPROP, which numerically computes the equilib¬ 
rium spectra of CR leptons and nuclei for given source 
functions and boundary conditions [9, 10]. I define a ref¬ 
erence model as follows. The source spectra are taken 
as broken power-law functions, qj oc {TZ/TZs)~'', with in¬ 
dex vi = 1.9 {v 2 = 2.38) below (above) TZb =9GV. The 
diffusion coefficient is taken as D{'R) = I3Dq {JZ/'Rq)^, 
with £>0 = 5 • 10^®cm^s“^, S =0.38, and TZo =4GV. 
The Alfven speed is va =33kms“^. The cylindrical dif¬ 
fusion region has radius r^nax = 30 kpc and half-height 
L = 3.9 kpc. A large nuclear reaction chain is set up, de¬ 
scribing the production of secondary j-type nuclei from 
fragmentation of heavier fc-type nuclei. The fragmen- 
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FIG. 1. (Color online) Elemental ratios B/C (top) and Be/B (bot¬ 
tom) from the reference model in comparison with the data [15- 
20]. The yellow bands are the estimated nuclear uncertainties (from 
Fig. 2). The red bands reflect the estimated parameter uncertain¬ 
ties for anticipated AMS data (from Fig. 3). 


tation rate is = l3kcnal^j{E), where rii are the 

number densities of the ISM nuclei (nn — 0.9 cin“^ and 
riHe =0.1cm“^) and is the production XS off i- 

type target at energy E. Under GALPROP, production 
XS's are evaluated from interpolation/fits to the data or 
from nuclear codes, such as CEM2k or LAQGSM, eventually 
normalized to the data [11-14]. 

Nuclear uncertainties — Data on isotopically sepa¬ 
rated targets and fragments have been collected by sev¬ 
eral experiments, but only in narrow energy ranges. A 
compilation is shown in Fig. 2 for the production of 

^Be, ®Be, and ^°Be isotopes from fragmentation of 
12 14,15 N and off hydrogen target at energy be¬ 
tween 30MeV/n and lOGeV/n. For Be production, I 
have also considered tertiary reactions such as B—^Be, 
j.e., those reactions where the progenitor nuclei are of 
secondary origin. All these processes account for > 90% 
of the Be and B production. At energy above than a 
few GeV/n, all the XS's are nearly constant in energy. 
The data in Fig. 2 are compared with XS parameteri- 
zations from GALPROP and from the popular formulae 
WNEW [32-34] and YIELDX [35]. Following earlier stud¬ 
ies [36, 37], the XS uncertainties are determiend by a 
re-fit of the GALPROP parameterizations (Tg{E) to the 
data. For each reaction, the corresponding XS are fit 
with the function cth(A) = aaG{bE), where a and b are 
free parameters representing normalization and energy 
scale. This procedure allows to determine a new set of 
XS's and their associated uncertainties corresponding to 


one-sigma confidence intervals. The uncertainty bands 
are shown in Fig. 2. The new XS's are often close to the 
original ac values, but the Be production in GALPROP 
is found to be over-estimated by a few percent. Such 
a Be overproduction was also reported in Aguilar et al. 
[15]. The estimated XS errors have been converted into 
uncertainties of the secondary/tertiary production terms 
and then propagated at Earth. Typical uncertainties are 
found to be ^ 5 % for B production and ^ 7-10% for Be 
production, with ^10% for ^°Be productions. The cor¬ 
responding uncertainties in the B/C and Be/B ratios are 
shown in Fig. 1 for the reference model (yellow bands). 

Modeling the AMS performance — I consider the B/C 
ratio at 2 - 200GeV/n and the Be/B ratio at 1 - 
lOOGeV/n. For these ratios, I compute the anticipated 
AMS data under the reference model. The number of 
j-type particles detected by AMS in each energy bin is 
estimated via the convolution of the CR flux with the 
detector acceptance [38], ANj = f ^jOjTjdE, where 
is the input spectrum, Qj is the total detector accep¬ 
tance and 7/ is the effective exposure time for a to¬ 
tal data taking period Tq. All input spectra are solar- 
modulated under the force-field approximation [39], us¬ 
ing (p = 550 MV for the AMS observation period. I con¬ 
sider the case of 10 bins per decade, log-uniformly spaced 
in energy, and a total exposure of QTq = 100 m^ sr day. 
The effective exposure time must account for the geo¬ 
magnetic field modulation which suppresses the Galactic 
GR flux below the cut-off rigidity, TZc ~ 0.5-20 GV, de¬ 
pending on the detector location. I adopt the Stormer 
model, TZc{t) = 20 GVp“^(t)cos^0M(t) [40], where 0 m( 0 
is the geomagnetic latitude and p{t) is the distance be¬ 
tween AMS and the geomagnetic dipole center in units 
of Earth’s radii. Their evolution depends on the ISS or¬ 
bit around the Earth. The function 7/ is computed as 
Tj{E) = Jj,^a{t)'Hj{t,E)dt, where a ~95% is the de¬ 
tector live-time, and 'Hj{t,E) is the geomagnetic trans¬ 
mission function, which is modeled as an 77-dependent 
smoothed step function, 7^ = [l -|- (77/77c)“^^] Its 
particle-dependence arises from the conversion TZ ^ E, 
while its time-dependence is contained in 77c (7). The 
integral Tj{E) has been numerically computed for all rel¬ 
evant isotopes by simulating 23,000 ISS orbits with pe¬ 
riod Tiss = 91 min and inclination 0iss =51.7°. Sys¬ 
tematic errors are assigned to be 1.5% for the B/C ratio 
and 1 % on the Be/B ratio, constant in the considered 
energy range [7]. From the estimated counts, the statis¬ 
tical errors associated with the B/C ratio are given by 
1/-\/A7 Vb + l/\/ANc, and similarly for the Be/B ratio. 


III. RESULTS AND DISCUSSIONS 

The AMS physics potential — The AMS capabilities 
in constraining the model parameters are first estimated 
without accounting for nuclear uncertainties. For this 
purpose, I have performed a scan in the parameter space 
Dq X L X va i>y running GALPROP 3,420 times over a 
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FIG. 2. (Color online) Fragmentation XS’s for ^Be, ^Be, and ^^Be production from B-C-N-0 collisions off hydrogen. The data 

are from A: Read & Viola 1984 [21], •: Webber et al. 1998 [22], □: Webber et al. 1990 [23], ▼: Olson et al. 1983 [24], ★: Fontes 1977 
[25], ■: Korejwo 2000 [26], 0: Korejwo et al. 2001 [27], A: Radin et al. 1979 [28], O: Ramaty et al. 1997 [29], X: Webber et al. 1998 
[30], Raisbeckl971 & Yiou 1971 [31]. The lines are from the WNEW (short-dashed), YIELDX (long-dashed), GALPROP (dotted), and 
the XS’s determined in this work (thick solid lines) with their uncertainty band. 
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FIG. 3. (Color online) Top: estimation of the AMS capabilities in constraining the parameters Dq, L, and va with with the B/C and 
Be/B ratios. Bottom: same as above after accounting for nuclear uncertainties in the Be-B production rates. 


19 X 15 X 12 grid. The resulting spectra are then tri- 
linearly interpolated to a 4x finer parameter grid, corre¬ 
sponding to 187,245 models. Hence, the B/C ratio pre¬ 
dicted by each model is tested against the artificial AMS 
data (generated with the reference model of Fig. 1), us¬ 
ing the method. The same procedure is done for Be/B 
ratio and for the two combined ratios. The results of are 
shown in Fig. 3, top panels, where the one-sigma contour 
regions are plotted as 2D projections of the parameter 
space. These regions are obtained from the surfaces 
of the B/C ratio and of the B/C-I-Be/B combination. The 
best-fit model on each plot (marked as “x”) always re¬ 
covers the true reference model. The complementarity of 
the two ratios in breaking the L-Dq degeneracy is appar¬ 
ent. While the B/C ratio constrains the two parameters 
into a tight region of the (L, Dq) plane, only the combi¬ 
nation B/C-t-Be/B allows to resolve their single values. 
The Alfven speed va is well determined by means of 
B/C data only. Tighter constraints be obtained using 
data below 2 GeV/n, provided that the solar modulation 
effect is well modeled. The accuracy of the measured 
parameters is 6Dq ~ 0.5 • 10^®cm^s“^, 6L ~ 0.5 kpc, 
and Sva ^ 2kin/s. This level of accuracy, from the esti¬ 
mated AMS capability, would represent quite a signihcant 
progress in CR propagation. 

Impact of nuclear uncertainties — The models con¬ 
strained by AMS are shown in Fig. 2 for both ratios (red 
bands). As clear from the figure, nuclear uncertainties 
(yellow bands) are dominating. In order to evaluate how 
these uncertainties affect the parameter reconstruction, 


I have repeated the parameter determination procedure 
after accounting for the XS errors in the x^ calculations. 
The results are shown in the bottom panels of Fig. 3. In 
comparison with the top panels, one can see that nuclear 
uncertainties have a dramatic impact on the parameters 
Dq and L. As shown in the hgure, the Dq/L degeneracy 
remain essentially unresolved when the nuclear uncer¬ 
tainties are taken into account. In fact, the information 
needed to break the Dq jL degeneracy is contained in the 
^°Be—^-^^B decay which produces only small variations in 
the Be/B ratio. Along with large uncertainties on the 
^°Be production, this information is also washed out by 
uncertainties in the more abundant ^’®Be and ^^B com¬ 
ponents of the Be/B ratio. At this point one may argue 
that a direct, ideal measurement of ^°Be at ~ 1-10 GeV/n 
would bring tighter constraints. Thus, I have repeated 
the calculations after considering XS uncertainties for the 
^°Be production only, be., assuming ideal knowledge of 
the other isotopes and infinite precision measurements. 
The sole uncertainties in the ^°Be production would limit 
the parameter reconstruction to 6Dq ^ 1.5 • 10^®cm^s“^ 
and 5L 1.5 kpc. This still represents a poor parame¬ 
ter determination in comparisons to the AMS potential. 
Nonetheless, given the current level of nuclear uncertain¬ 
ties, a direct measurement of ^°Be flux (even if affected 
a few % systematic errors) would probably bring better 
information than a precise Be/B measurement. 

Single-reaetion XS bias — It is instructive to study 
the dependence of the best-fit parameters on single XS 
reactions. An example is the anti-correlation between 
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FIG. 4. (Color online) Best-fit parameters Dq? L, and Va and the 
corresponding Xmin ^ function of the relative bias introduced in 
the XS normalization for the reactions and 


the reaction and the parameter L [41]. Here 

I consider the reactions ^-^^B and 7>^Be. After 
the introduction of systematic biases in their XS normal¬ 
izations, I have generated new reference-mo del predic¬ 
tions and new AMS mock data. Then, I have repeated 
the parameters reconstruction procedure using the nom¬ 
inal set of unbiased XS’s. The results are shown is in 
Fig. 4. The best-fit parameters are plotted as a function 
of the bias induced on the XS normalization for the con¬ 
sidered reactions. The horizontal dotted lines show the 
true parameter values. In this case, the best-fit param¬ 
eters are not correctly reconstructed. While the deter¬ 
mination of va appears rather stable, the parameter L 
and Dq are mis-inferred by large factors after the intro¬ 
duction of a few % deviation in the relevant XSs. This 
shows, again, that the constraining power of the Be/B 
ratio is weakened by uncertainties in the ^^B or ^’®Be 
production, rather than by the large errors on the ^°Be 
production. It is also interesting to look at the evolution 
of the best-fit Xmin function of the bias. This is shown 
in the bottom panels of the figure. The best of the 
B/C ratio appears insensitive to XS biases. In fact the 
B/C ratio is approximately given by B/C oc rB/{D/L), 


i.e., any deviation in the B production rate Tg can be 
re-absorbed by a different determination of Dq/L. This 
is not the case for the Be/B ratio because its high-energy 
plateau (at > 10 GeV/n) is almost independent on prop¬ 
agation effects. Hence, discrepancies between Be/B data 
and model predictions cannot be re-absorbed by the fit in 
term of different parameter combinations: they can only 
arise from nuclear physics inputs. An example of this 
is found in Aguilar et al. [15], where the small discrep¬ 
ancy between the Be/B data and the model predictions 
was ascribed to the XS’s for Be production. The Be/B 
ratio, like other secondary-to-secondary ratios [42], may 
be therefore used as a diagnostic tool to detect possible 
biases in the production XS’s. 


IV. CONCLUSIONS 

My estimates show that the AMS experiment can 
provide tight constraints on the key parameters Dq, L, 
and VA- Given the level of precision expected by AMS, 
nuclear uncertainties in secondary production models 
are found to be a major limitation in the interpretation 
of secondary GR nuclei. Once nuclear uncertainties 
are accounted, the Dq/L degeneracy remains poorly 
resolved. With the current status of nuclear data, the 
Be/B ratio appears not to bring valuable information 
for the parameter extraction. Isotopically resolved CR 
measurements such as the ^°Be/®Be ratio are prefer¬ 
able, though the ^°Be production rate is affected by 
large uncertainties. On the other hand, precise data 
on the Be/B ratio at E >I0GeV/n may represent a 
powerful tool to test the nuclear physics inputs of the 
propagation models, and in particular to detect possible 
biases in single reactions that may cause a parameter 
mis-determination. It also worth stressing that this 
problem has a direct impact in dark matter searches. 
My study provides a concrete case study for Gondolo’s 
plea to the nuclear physics community [43] . In summary, 
nuclear uncertainties are a major limiting factor for 
further progress in GR propagation. The collection 
of new nuclear data, within a dedicated program of 
XS measurements and modeling, would enable to fully 
exploit the potential of the AMS data. 
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